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Past studies 1 3 suggest that long-duration 7-ray bursts (GRBs) have a stan- 
dard energy of E 1 ~ 10 51 erg in ultra-relativistic ejecta when corrected for 
asymmetry ("jets"). However, recently 2,3 a group of sub-energetic bursts, in- 
cluding the peculiar GRB 980425 associated 4 with SN 1998bw (E 1 « 10 48 erg), 
has been identified. Here we report radio observations of GRB 030329, the 
nearest burst to date, which allow us to undertake calorimetry of the ex- 
plosion. Our observations require a two-component explosion: a narrow (5°) 
ultra-relativistic component responsible for the 7-rays and early afterglow, 
and a wide, mildly relativistic component responsible for the radio and opti- 
cal afterglow beyond 1.5 days. While the 7-rays are energetically minor, the 
total energy release, dominated by the wide component, is similar 1 3,5 to that 
of other GRBs. Given the firm link 6 7 of GRB 030329 with SN 2003dh our 
result suggests a common origin for cosmic explosions in which, for reasons 
not understood, the energy in the highest velocity ejecta is highly variable. 

We initiated observations of the nearby GRB 030329 (z = 0.1685) in the centimetre 
band approximately 13.8 hours after the burst. The log of the observations and the 
resulting lightcurves are displayed in Tables 1 and 2 and Figure 1. The afterglow was also 
observed extensively in the millimetre (100 GHz) and sub-millimetre (250 GHz) bands. 8 
While this is the brightest radio afterglow detected to date, the low redshift results in a 
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peak luminosity, L UyP (8.5 GHz) « 1.8 x 10 31 erg s -1 Hz -1 , typical 9 of other long-duration 
GRBs. 

The observed rapid decline, F u oc t -19 at t ^ 10 d and the decrease in peak flux at v ^ 
22.5 GHz (Figure 1) are the hallmarks of a collimated explosion. In this framework, 10 the 
sharp decline (or "jet break") occurs at the time, tj, when T(tj) ~ Qj 1 due to relativistic 
abberation ( "beaming" ) and rapid side-ways expansion; here Y is the bulk Lorentz factor 
and Qj is the opening angle of the jet. We model the afterglow emission (cf. ref. 5,11) 
from 4.9 to 250 GHz assuming a uniform 10 as well as a "wind" 12 (particle density profile, 
p oc r~ 2 , where r is the distance from the source) circumburst medium. Neither model is 
strongly preferred, but t, ira d ~ 9.8 d is required (Figure 1). 

Using the inferred particle density of n ~ 1.8 cm -3 and assuming a 7-ray efficiency, 
e 7 = 0.2 (see ref. 3) we infer #, irad ~ 0.3 rad, or 17°. The kinetic energy in the explosion 
corrected for collimation is E K = fbE K:iso « 2.5 x 10 50 erg, where /& = [1 — cos(^)] is 
the beaming fraction and E Kjiso is the isotropic equivalent kinetic energy. This value is 
comparable to that inferred from modeling of other afterglows. 5 

In contrast to the above discussion, Price et al. 13 note a sharp break in the optical 
afterglow at t = 0.55 d (Figure 2). The X-ray flux 14 tracks the optical afterglow for the 
first day, with a break consistent with that seen in the optical. Thus the break at 0.55 d 
is not due to a change in the ambient density since for typical parameters 15 ' 16 the X-ray 
emission is not sensitive to density. However, unlike the optical emission the X-ray flux 
at later times continues to decrease monotonically. Thus we conclude that there are two 
emitting components: one responsible for the early optical and X-ray emission and the 
other responsible for the optical emission beyond 1.5 days. 

The first component, given the characteristic t~ 2 decay for both the X-ray and optical 
emission, is reasonably modeled by a jet. For the parameters used above (n, e 7 ) the 
opening angle is 0.09 rad or 5°. 

The resurgence in the optical emission at 1.5 d requires a second component. An 
increase in the ambient density cannot explain this resurgence since the predicted de- 
crease in radio luminosity, arising from the increase in synchrotron self-absorption, is not 
observed (Figure 1). An increase in the energy of the first component, for example by 
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successive shells with lower Lorentz factors as advocated by Granot et al., 17 is ruled out 
by the lack 8 of strong radio or millimetric emission expected 18 from reverse shocks. 

Thus, by a process of elimination, we are led to a two-component explosion model 
in which the first component (a narrow jet, 5°) with initially larger T is responsible for 
the 7-ray burst and the early optical and X-ray afterglow including the break at 0.55 
d, while the second component (a wider jet, 17°) powers the radio afterglow and late 
optical emission (Figure 2). The break due to the second component is readily seen in 
the radio afterglow, but is masked by SN2003dh in the optical bands, thus requiring 
careful subtraction (Figure 2). Such a two-component jet finds a natural explanation in 
the collapsar model. 19 

The beaming-corrected 7-ray energy, emitted by the narrow jet, is only E 1 R3 5 x 10 49 
erg, significantly lower than the strong clustering 3 around 1.3 x 10 51 erg seen in most 
bursts. Similarly, the beaming-corrected X-ray luminosity 14 at t — 10 hours, a proxy 
for the kinetic energy of the afterglow on that timescale, is Lx,io ~ 3 x 10 43 erg s~ 4 , 
a factor of ten below the tightly clustered values 2 for most other bursts. However, the 
second component, which is mildly relativistic (as determined by the lower energy peak 
of its spectrum), carries the bulk of the energy, as indicated by our modeling of the radio 
emission. We note that our model, with the energy in the lower Lorentz factor component 
dominating over the narrow ultra- relativistic component, is not consistent with "universal 
standard jet" model. 20 

The afterglow calorimetry presented here has important ramifications for our under- 
standing of GRB engines. Recently, we have come to recognize a sub-class of cosmological 
GRBs marked by rapidly fading afterglows at early time (i.e. similar to GRB 030329). 
These events are sub-energetic 2 ' 3 in E y and early X-ray afterglow luminosity. However, 
as demonstrated by our calorimetry of GRB 030329, such bursts may have total explosive 
yields similar to other GRBs (Figure 3). 

This leads to the following conclusions. First, radio calorimetry, which is sensitive to 
all ejecta with T ^ few, shows that the explosive yield of the nearest "classical" event, 
GRB 030329, is dominated by mildly relativistic ejecta. Ultra-relativistic ejecta which 
produced the 7-ray emission is energetically unimportant. Second, the total energy yield of 
GRB 030329 is similar to those estimated for other bursts. Along these lines, the enigmatic 
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GRB 980425 associated 4 with the nearby supernova SN 1998bw also has negligible 7-ray 
emission, -E 7)iso ~ 8 x 10 47 erg; however, radio calorimetry 21 shows that even this extreme 
event had a similar explosive energy yield (Figure 3). The newly recognized class of 
cosmic explosions, the X-ray Flashes, 22 exhibit little or no 7-ray emission but appear to 
have comparable X-ray and radio afterglows to those of GRBs. Thus, the commonality 
of the total energy yield indicates a common origin, but apparantly the ultra-relativistic 
output is highly variable. Unraveling what physical parameter is responsible for the 
variation in the "purity" (ultra-relativistic output) of the engine appears to be the next 
frontier in the field of cosmic explosions. 
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6.08 ± 0.10 


4.48 ±0.09 


2.82 ±0.21 


2.84 ±0.20 




June 4.01 


66.53 


1.94 ±0.06 


6.20 ± 0.08 


4.93 ±0.06 




2.56 ±0.12 





Table 1. Radio observations made with the Very Large Array (VLA) and the Australia Tele- 
scope Compact Array (ATCA). Observations commenced on March 30.06 UT with a single 
7-hour observation with ATCA on Mar. 30.53 UT. In the initial observation we detected a point 
source at a(J2000)=10 h 44 m 49.95 s , <5(J2000)=21°31'17.38", with an uncertainty of about 0.1 
arcsec in each coordinate, consistent with the position of the optical counterpart. In all VLA 
observations we used the standard continuum mode with 2 x 50 MHz bands. At 22.5 and 43.3 
GHz we used referenced pointing scans to correct for the systematic 10 — 20 arcsec pointing 
errors of the VLA antennas. We used the extra-galactic sources 3C 147 (J0542+498) and 3C 286 
(J1331+305) for flux calibration, while the phase was monitored using Jllll+199 at 1.43 GHz 
and J1051+213 at all other frequencies. The ATCA observations were performed at 4.80, 6.21, 
8.26, and 9.02 GHz with a bandwidth of 64 MHz in each frequency. The phase was monitored 
using J1049+215, while the flux was calibrated using J1934-638. The data were reduced and 
analyzed using the Astronomical Image Processing System (VLA) and the Multichannel Image 
Reconstruction, Image Analysis and Display package (ATCA). The flux density and uncertainty 
were measured from the resulting maps by fitting a Gaussian model to the afterglow. In addition 
to the rms noise in each measurement we estimate a systematic uncertainty of about 2% due to 
uncertainty in the absolute flux calibration. 
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Apr 1.12 


2.64 


18.29 ±0.28 
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13.60 ±0.65 
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16.75 ±0.27 
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11.78 ±0.52 


Apr 3.07 
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20.36 ±0.45 
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10.35 ±0.49 
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36.32 


9.15 ±0.50 
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28.69 ±0.26 


May 6.83 


38.35 


7.87 ±0.50 
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8.58 


29.29 ±0.49 


May 8.73 


40.25 


6.70 ±0.50 


Apr 7.89 


9.41 


29.15 ±0.44 


May 10.76 


42.28 


6.49 ±0.50 


Apr 9.89 


11.41 


30.78 ±0.51 


May 15.76 


47.28 


5.74 ±0.50 


Apr 11.05 


12.57 


28.52 ±0.51 


May 20.70 


52.22 


5.69 ±0.53 


Apr 11.88 


13.40 


29.92 ± 0.44 


May 22.76 


54.28 


4.78 ±0.78 


Apr 13.05 


14.57 


27.90 ± 0.44 


May 24.76 


56.28 


4.31 ±0.55 


Apr 13.87 
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24.74 ± 0.44 


May 25.56 
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5.04 ±0.84 
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23.60 ±0.32 


May 26.75 
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3.99 ±0.63 


Apr 16.96 
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23.06 ±0.24 
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3.96 ±0.58 
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20.51 ±0.24 


May 29.82 


61.34 


4.35 ±0.50 


Apr 19.95 


21.47 


19.27 ±0.38 


May 30.76 


62.28 


2.65 ±0.72 


Apr 20.72 


22.24 


17.53 ±0.33 


June 2.54 


64.06 


3.13 ±0.76 



Table 2. Radio observations at 15.3 GHz made with the Ryle Telescope at Cambridge (UK). 
All observations were made by interleaving 15 minutes scans of GRB 030329 with 2.5 minutes 
scans of the phase calibrator J1051+2119. The absolute flux scale was calibrated using 3C 48 and 
3C 286. We used 5 antennas providing 10 baselines in the range 35 - 140 m. Since the position 
of the source is well known the in-phase component of the vector sum of the 10 baselines was 
used as an unbiased estimate of the flux density. The typical rms fluctuation on the signal in a 
32-s integration period is approximately 6 mjy. We also add a systematic uncertainty of about 
2% due to uncertainty in the absolute flux calibration. 
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Figure 1. Radio lightcurves of the afterglow of GRB 030329. All measurements include la error 
bars which in most cases are smaller than the symbols. The data are summarized in Tables 1 
and 2. The solid lines are models of synchrotron emission from collimated relativistic ejecta 
expanding into uniform 10 (thick) and wind 12 p oc r~ 2 (thin) circumburst media; these models 
include the millimetre and sub-millimetre data. 8 We find Xr = 31-3 and 39.8 (164 degrees of 
freedom) for the uniform density and wind models, respectively; these include a 2% systematic 
error added in quadrature to each measurement. The large values of Xr are dominated by 
interstellar scintillation (ISS) at v <^ 15 GHz and mild deviations from the expected smooth 
behavior at the high frequencies. Comparing the data and models, we find rms flux modulations 
of 0.25 at 4.9 GHz, 0.15 at 8.5 GHz, and 0.08 at 15 GHz, as well as a drop by a factor of three 
in the level of modulation from ~ 3 to 40 days. These properties are expected in weak ISS as 
the fireball expands on the sky. The inferred source size of about 20 ^as (i.e. ~ 2 x 10 17 cm) 
at t ~ 15 days is in close agreement with theoretical expectations. 23 In the uniform density 
model the jet break occurs at t ~ 10 d corresponding to an opening angle, 9j ~ 0.3 (17°). 
From the derived synchrotron parameters (at t = tj): v a ~ 19 GHz, v m ~ 43 GHz, F U)Tn ~ 96 
mJy we find an isotropic kinetic energy, E K ^ so m 5.6 x 10 51 

^ 'c 13 er §) a circumburst density 
n = 1.8^^3 cm -3 , and the fractions of energy in the relativistic electrons and magnetic field of 
0.16^^3 and O.lOz-^^ 4 , respectively; here v c = 10 13 ^ Cj i3 is the synchrotron cooling frequency, 
and a constraint on Inverse Compton cooling as advocated by Sari & Esin 24 indicates is c ,i3 ^ 1- 
The beaming-corrected kinetic energy is Ek ~ 2.5 x lO 50 ^^ erg, typical of other well-studied 
long-duration GRBs. 5 The parameters derived from the wind model are consistent with those 
from the uniform density model to within 10%. 
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10" 1 10° 10 1 10 2 

Time Since the Burst (days) 

Figure 2. Radio to X-ray lightcurves of the afterglow of GRB 030329. The optical data, 
from Price et al. 13 and the GRB Coordinates Network, 25-27 have been corrected for Galactic 
extinction, Ar = 0.067 mag; we note that the latter are preliminary data. The dotted line 
is the model proposed by Price et al. 13 for the early optical emission ? with t^opt « 0.55 d. 
The dashed line is an extapolation of the uniform density model presented in Figure 1 to the 
optical -R-band with 1^0,13 = 2; this value is somewhat larger than the rough limit discussed in 
Figure 1 but may be consistent with the uncertainty in the model parameters. The model in 
the X-ray band is based on the measured 14 optical to X-ray spectral slope and an extapolation 
of the uniform density model presented in Figure 1. The sharp increase in the optical flux at 
t <^ 1.5 d is due to the deceleration of the slower second jet component. Finally, the dot-dashed 
line is the optical emission from SN1998bw at the redshift of GRB 030329, z = 0.1685, used 
as a proxy for SN 2003dh. 6 The solid line in the i?-band is a combination of the SN and the 
two jet components, whereas in the radio and millimetre bands it is the uniform density model 
presented in Figure 1. In the X-ray band the model is dominated by the narrow jet component. 
While this two-component jet model provides a reasonable fit to the data, there are still some 
discrepancies which can be resolved by accurate photometry and a more careful subtraction of 
SN 2003dh. The latter will also allow a more precise determination of the putative second jet 
break in the optical band at i^rad- 
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Figure 3. Histograms of various energies measured for GRBs. We plot the beaming-corrected 
7-ray energy, 3 E^, the kinetic energy inferred from X-rays at t = 10 hr, 2 Ek,x, and total 
relativistic energy, E rc \ = E 7 +Ek, where Ek is the beaming-corrected kinetic energy inferred 21,5 
from the broad-band afterglow. The energy in X-rays, Ek,x = Lxt/e e {ax — 1), with t = 10 
hr, e e = 0.1, and ax = 1-3 is the median decay rate in the X-ray band. For GRB 980519 we 
find that the evolution of the radio emission requires a much wider jet, 9j ~ 0.3, than what is 
inferred from the optical, 9j ~ 0.05; here we assume z = 1. We therefore infer Ek ~ 2 x 10 50 
erg from the radio data compared to E^ ~4x 10 erg. The 7-ray energy of GRB 980425 is an 
upper limit since the degree of collimation is not known. For the kinetic energy we use the value 
derived by Li & Chevalier 21 based on the radio evolution of SN 1998bw. There is a significantly 
wider dispersion in £L and Ek,x as compared to the total explosive yield. This indicates that 
engines in cosmic explosions produce approximately the same quantity of energy thus pointing 
to a common origin, but the quality of these engines, as indicated by ultra-relativistic output, 
varies widely. 



